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Deterministic compressed sensing matrices: 
Construction via Euler Squares and applications 

R. Ramu Naidu, Phanindra Jampana and C. S. Sastry 


Abstract —In Compressed Sensing the matrices that satisfy the 
Restricted Isometry Property (RIP) play an important role. But 
to date, very few results for designing such matrices are available. 
For applications such as multiplier-less data compression, binary 
sensing matrices are of interest. The present work constructs 
deterministic and binary sensing matrices using Euler Squares. 
In particular, given a positive integer m different from for a 
prime p, we show that it is possible to construct a binary sensing 
matrix of size m x c(mfi)^, where p is the coherence parameter 
of the matrix and c € [1,2). The matrices that we construct 
have small density (that is, percentage of nonzero entries in the 
matrix is small) with no function evaluation in their construction, 
which support algorithms with low computational complexity. 
Through experimental work, we show that our binary sensing 
matrices can be used for such applications as content based 
image retrieval. Our simulation results demonstrate that the 
Euler Square based CS matrices give better performance than 
their Gaussian counterparts. 

Index Terms —Compressed Sensing, Coherence, RIP, Binary 
sensing matrices, Euler Squares, CBIR. 


I. Introduction 

Compressed Sensing (CS) aims at recovering high dimen¬ 
sional sparse vectors from considerably fewer linear mea¬ 
surements. The problem of sparse recovery through Iq norm 
minimization is not tractable. E. Candes JTOl . D. Donoho IITtI . 
have made pioneering contributions reposing the problem as a 
simple Linear Programming Problem (LPP). They have then 
established the conditions that ensure the stated equivalence 
between the original Iq problem and its reposed version. It is 
known that Restricted Isometry Property (RIP) is one sufficient 
condition to ensure the equivalence. Random matrices with 
Gaussian or Bernoulli entries have been shown to satisfy RIP 
with high probability Q. 

In the recent literature on CS 0, HI, 0, the deterministic 
construction of CS matrices has gained momentum. A. Amini 
et. al. im, R. Devore Ill3l . P. Indyk im and S. Li et. al. 
ED have constructed deterministic binary sensing matrices 
using ideas from algebra and graph theory. In the present work, 
however, we attempt to construct deterministic binary sensing 
matrices using Euler Squares. The advantages of the proposed 
methodology are as follows: 
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• Matrices of general row size (different from a prime and 
its square) and small density are constructed. 

• Simplicity in construction is achieved 

Here, density is defined as the ratio of number of nonzero 
entries to the total number of entries of the matrix. Sparse 
sensing matrices may contribute to fast processing with low 
computational complexity in compressed sensing ina. 

In the recent past, Gaussian and ±1 Bernoulli random 
matrices (which are shown to satisfy RIP ID) have been used 
to project data into lower dimension for the purpose of clas¬ 
sification ll26l . There are, however, the following advantages 
of using deterministic binary matrices for the stated purpose: 

• Binary matrices being sparse and possessing 0,1 as 
elements provide multiplier-less and faster dimensionality 
reduction operation, which is not possible with Gaussian 
matrices 

• There is a nonzero probability of non-compliance with 
RIP for Gaussian matrices 

The problem (ED, l28l ) of searching for similar images in 
a large image repository based on content is called Content 
Based Image Retrieval (CBIR). CBIR has several multimedia 
applications like text based retrieval (such as google search) 
and retrieval from medical databases. As images of large 
size typically involve more complexity, one needs to project 
them to lower dimensional spaces. It is demonstrated that 
the proposed binary sensing matrices project data into lower 
dimensional spaces in such a way that the reduced vectors 
are useful for the purpose of CBIR. Moreover, the proposed 
dimensionality reduction technique through binary sensing 
matrices allows for reconstruction, which is important for 
tele-medicine (ID, ED). The other dimensionality reduction 
techniques and even the sparsity seeking Dictionary based 
methods El in general do not provide this advantage. 

The paper is organized into several sections. In section [III 
basic CS theory and the equivalence between Iq norm and 
li norm problems are given. In section [Bll the deterministic 
construction procedure of binary sensing matrices using Euler 
Squares is presented. Section HVl gives comparison with exist¬ 
ing constructions. In sectionlV] an application to content based 
image retrieval is demonstrated. Concluding remarks are given 
in section m 

II. Sparse recovery from linear measurements 

As stated already, CS refers to the problem of reconstruction 
of an unknown vector u G from linear measurements 
y = ((m, (^i),..., (u, 0m)) C with (m, 0^) being the inner- 
product between u and 0j. The basic objective in CS is to 
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design a recovery procedure based on the sparsity assumption 
on u when the number of measurements m is very small 
compared to M. Sparse representations have merit for various 
applications llT2l . ifTSl . lIZTil . If30l . ll^ in areas such as 
image/signal processing and numerical computation. 

A vector u G is /c—sparse if it has at most k nonzero 
coordinates. Let ||u||o stand for |{i | Vi ^ 0}|. The problem of 
obtaining the sparse vector from its linear measurements may 
be posed as 

Pq : min ||w||o subject to = y. 

V 

This Zo~miiiimization problem is NP-hard Cl in general. 
There are greedy algorithms for solving Pq problem and 
orthogonal matching pursuit (OMP) is one of the popular 
methods C3. Several researchers (for example, D. Donoho 
d and E. Candes ifTOl ) have established the conditions that 
ensure the recovery of solution to Pq from : 

Pi : min ||u||i subject to fhu = y. 

V 

Here ||u||i denotes the (i-norm of the vector v G Denote 
the solution to Pi by /$(y) and solution to Pq by u^(y) G 


A. On the equivalence between Pq and Pi problems 


The coherence /i($) of a given matrix $ is the largest 
absolute inner-product between different normalized columns 
of $. Denoting the fc-th column in $ by (pk, one defines 

the coherence, as /i($) = ^ 


max 


1< i,j< M, II 2 


. For a 


M-r. 


called the 


matrix of size m x M, p, satisfies: p > 

Welch bound. It is known ca that for /^-coherent matrices 
<1), one has u%{y) = f^{y) = u, provided u is fc—sparse with 

fc<i(l + i). 


The Restricted Isometry Property (RIP) plays an important 
role lITOl in CS as it establishes the equivalence between the Pq 
and Pi problems. An mxM matrix $ (which we refer to a CS 
matrix) is said to satisfy the Restricted Isometry Property(RIP) 
of order k with constant 6k if for all fc—sparse vectors x G 
we have 


{l-Sk)\\x\\l<\\^x\\l<il + Sk)\\x\\l. (1) 

It is known caci that the RIP along with some conditions 
on Sk imply the equivalence between the Pq and Pi problems. 

One of the important problems in CS theory deals with 
constructing CS matrices that satisfy RIP for the largest 
possible range of k. It is known that the random constructions 
satisfy the RIP for the largest possible range on sparsity, which 
= 0 . 


B. Existing deterministic constructions 

To the best of our knowledge, designing good determin¬ 
istic constructions of RIP matrices is still a very interesting 
problem. R. Devore lfT3l has constructed deterministic binary 
sensing matrices of size x with coherence where 
p is a prime power and 0 < r < p. Later on, S. Li, F. Gao 
et. al. 12111 have generalized the work in ca, constructing 


binary sensing matrices of size \V\q x where q is any 

prime power and V is the set of all rational points on algebraic 
curve X over finite field F^. P. Indyk ifTSll has constructed 
binary sensing matrices using hash functions and extractor 
graphs with sizes iogn)°(i) X n, where r <C n. A. 

Amini et. al. 0 have constructed binary sensing matrices 
using OOC codes. In all these constructions, CS matrices are 
given for specific sets of row sizes. J. Bourgain et. al. 13 
have constructed RIP matrices of order k > 7712 +'^, for some 
e > 0 and M^~'^ < m < M using additive combinatorics. 
It is remarkable that this construction overcomes the natural 
barrier k — 0{m2) for those based on coherence. J. L. 
Nelson et. al l29l have given lower and upper bounds on 
maximal possible column size in terms of fixed row size m 
and coherence p. They have given a class of deterministic CS 
matrices of size p x p^ with coherence 2 ^, where p is a 
prime number. R. Calderbank et.al. 0, 0 have constructed 
CS matrices of size 2* x 2^’'+^)^ with coherence 2’'“ 2 , using 
the Delsarte-Goethals codes, where I is an odd number and 
0 < r < is a constant integer. M. A. Herman et.al. 03 
have constructed Gabor frame of size p x p^ with coherence 
using Alltop sequence, where p > 5 is any prime number. 
The constructions, possessing sparsity k = y/m have been 
given by Calderbank et. al. in fS) based on tight frames such 
as TO X TO^ chirp matrices and All top Gabor frames 13, here 
TO is prime or prime power. The authors in 122, ll34ll have 
constructed matrices which achieve the square root bottleneck 
asymptotically. 

In the present work, binary sensing matrices are constructed 
using Euler Squares. In particular, given any positive integer 
TO ^ p,p^ (for prime p), the procedure allows to construct 
binary CS matrix of size to x M, where M is c{mp)^, with 
p being the coherence parameter of the matrix (which is 
dependent on to as explained in sub-sequent sections) and 
c G [1,2). As the proposed methodology does not involve any 
function evaluations and gives matrices of small density, this 
method involves less complexity as compared to some of the 
existing ideas stated above. The applicability of construction 
is demonstrated through an application to CBIR. 

III. Euler Squares for constructing CS matrices 
A. Euler Squares 

An Euler Square of order n, degree k and index n, fc is a 
square array of A:—ads of numbers, (a^i, 0 ^ 2 , ■ • •, flyfe), 
where G {0,1, 2,..., n — 1}; r = 1, 2,..., fc; i, j = 
1,2,..., 72, Ti ^ A:, CLipj‘ ^ cl'iqj- and Qjpjj’ ^ Ojqjy for p ^ q 
and {a^jr + + 1) ^ (opgr + l)(ap 9 s + 1) for 7 ^ p and 

i ^ 9 - 

Harris F. MacNeish 123 has constructed Euler Squares for 
the following cases: 

1. Index p,p — 1, where p is a prime number. 

2. Index p’’,p’^ — 1, for p prime. 

3. Index n,k, where n = 2^p]^p2^ ■ ■■ tPi' for distinct odd 
primes pi,p 2 ,... ,p/. Here, A: -f 1 equals the least of the 
numbers 2\p[\pl^,... ,p\’-. 

For example, the Euler Square of index 3, 2 is as follows: 
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0,0 

1,1 

2,2 

1,2 

2,0 

0,1 

2,1 

0,2 

1,0 


Lemma 1. ^2?l/ Let k' < k. Then the existence of the Euler 
Square of index n, k implies that the Euler Square of index 
n, k' exists. 


B. Deterministic constructions via Euler Squares 

Using the concept of Euler Square of index n, k, deterministic 
binary CS matrices $ that possess small coherence can be 
constructed. In particular, for n > 3,fc > 2, aCS matrix of 
size nkxnf, which may be treated as a block matrix consisting 
of k number of n x nf blocks can be obtained. Each column 
in <1) corresponds to a fc-ad in the Euler Square. The matrix $ 
is defined as follows: Eor 1 < i < nk, 1 < j < n^. 


4>ij 


1 if = i — l(mod n) 

0 otherwise, 


( 2 ) 


where (a^) is the A:—ad, (oj); is element in fc—ad 

and [a;J is the largest integer not greater than x. To summarize, 
for j = 1, 2,..., the column fj of $ is generated as 
an nk— binary vector from the k-ad (aj) with 1 occurring 
at the positions (Z — l)n + {{aj)i + 1) for Z = 1,2,..., fc. 

There are exactly fc ones in each column of $ and size of 
the matrix is nk x n^. Eor example, the 6x9 matrix given by 
the Euler Square of index 3, 2 is as follows: 


/ 1 0 0 0 0 

0 10 10 
0 0 10 1 
1 0 0 0 1 

0 10 0 0 

\ 0 0 1 1 0 


1 0 1 0 \ 

0 0 0 1 

0 10 0 

0 0 0 1 

110 0 
0 0 10/ 


The following lemma finds the bound on the coherence of $. 


Lemma 2. The coherence of $ is at most equal to p 

Proof: Suppose there exist two columns fc, fd of $ 
such that the overlap between them is at least 2. That is, if 
Uc,Ud are the support vectors of fcfd, then \ucr)Ud\ > 
2, which implies that there exist two fc—ads (oc) = 
{Oijl^CLij 2 , ■ • . ,CLijkf (ttd) — (UpQ, Up/ 2 , ■ • ■ , O.plk') with i = 
+1. nandp= +1, Z = d- [^\ n 

such that 

Uijr = Opic and a^s = Op/g for some 1 < r, s < fc. (3) 

Case 1: Suppose i = p and j Z, from Q, Op^y = Op/r 
and Oijs = aus for j I, which is a contradiction from the 
definition of Euler Square. 

Case 2: Suppose i p and j = 1. Erom ©, aur = a-pir and 
aijs = Upjs for i p, which is again a contradiction. 

Case 3: Suppose i p and j ^ Z, from (3), Oijr = Up/r and 
Uzjs = cLpis- Since, Oijr, cLpis € {0,1,..., 7t 1} for all 1 


r,s < k, we have a^y + 1 = Up/r + 1 and a^y + 1 = Op/s + 1, 
which implies that {aijr+l){aijs + f) = (ap/r.+l)(ap/s+l) for 
i p and j ^ Z, which is a contradiction from the definition 
of Euler Square. 

Hence there are no two such fc—ads satisfying ([3l. There¬ 
fore, the overlap between any two columns of $ is at most 1. 
Since each column in $ contains fixed number of fc ones, it 
follows that the coherence of $, /r($), is at most equal to p 

■ 

Remark 1: The maximum possible column size of any binary 
matrix is JU, where m is the row size, fc is the number 
of ones in each column and r — 1 is the overlap between any 
two columns. Euler Square of index n, fc results in a binary 
sensing matrix of size nk x n^. In this matrix each column 
contains fc number of ones and coherence is at most f. The 

k 

maximum possible column size is thus 

i^ = 0(n2f]=0((m/r)2). (4) 

121 

Eor the Euler Square based matrix m = nk,M = n^ and 
fi = J, and hence M = which is in the order of 

maximum possible column size. Hence the aspect ratio is also 
in the maximum possible order. The density for these matrices 

- I 

The following proposition Cl relates the RIP constant Sk' and 

p. 

Proposition 3. Suppose that fi,..., fm <dre the unit norm 
columns of the matrix $ possessing the coherence p,. Then <!> 
satisfies the RIP of order k' with constant Sk' = {k' — \)p. 

Erom lemma |2] and Proposition [3 it follows that the matrix $ 
so constructed satisfies RIP. 

Theorem 4. The matrix $o = satisfies RIP with Sk' = 
^ for any k' < k + 1. 

Remark 2: m = nk, M = n"^ and fc' < fc -I- 1 gives 
fc = ^,n = \/M which implies that fc = Consequently 
k'im^M) <^p+l. 

The previous procedure can be generalized for any row 
size m different from a prime and its square. The following 
theorem summarizes the main result. 

Theorem 5. Suppose m is any positive integer different from 
p, p^ for a prime p. Then there exists a binary sensing matrix 
of size m X M with coherence u = 

' m 

Proof: Case-1: If m = p*, i > 2 then m can be written 
as m = Since Z > 2, we have (ff~^ — 1) > p. Prom 

1241, it is known that the Euler Square of index — 1 

exists. Since — 1) > p, from lemma [T] Euler Square of 
index p®“^,p exists. If we construct the binary matrix using 
this Euler Square, then the row size of it becomes p*. 

Case-2: Let m be any integer such that m p'‘ for i > 2. 
Prom the fundamental theorem of arithmetic m is factorized as 
m = 2'^pf^plf ■ ■ - pV■ Let fc' = min{ 2 ’’,p^\p 2 ^,... ,pp} and 

*a = 0(b) implies that, there exist two constants ci,C 2 such that cib < 
a < C 2 b. 
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k = k' — 1. With out loss of generality assume that k' = 2’’ 
(The following arguments hold true even if k' = p''/ for some 
s € {1,2,..., (}). From the construction of MacNeish UtM . it 
is known that the Euler Square of index m, k exists. Let mi = 
PiP 2 ^ ■ ■ - Pi' ^nd without loss of generality assume that = 
min{p{\p 2 ^,... ,pp}. Since 2’" < (j>Y — 1), Euler Square of 
index mi, {p[^ — 1) exists, which in view of lemma [1] implies 
that the Euler Square of index mi, 2'’ exists. The binary CS 
matrix constructed through this Euler Square has a row size 
of mi2'’, which is m. ■ 

Remark 3; In both the cases above M = (mp,)^, which is 
in the order of maximum possible column size as stated in 
Remark-1 and they allow for successful signal recovery of 
sparsity k'{m,M) < + 1. 

The afore-described construction methodology results in a 
binary CS matrix of size nkxn^ from an Euler Square of index 
n, k. As explained in Remark 1, the maximum possible column 
size for our construction in terms of row size nk and coherence 

is = 0(n^). Since, < 2n^, we have ci = 1 

and C 2 = 2. In Theorem 5, we give the construction for c = 1 
and in the next section, when n is of the form 2 ’'p{^p 2 ^ ■ ■ - pY 
(that is, n is not a power of single prime) we try to increase the 
column size of above construction to c(n^), where c G (1, 2), 
leaving the row size and the coherence intact.The procedure 
discussed above uses smallest factor, while in increasing the 
column size, other factors are also used. 

C. Column extension of CS matrices constructed from Euler 
Squares 

As before, let n = 2^pf p'f-..., p^f, k'= 

min{2’',p{\p2^,... ,pp} and k = k' — 1. Suppose 

is the matrix of size nk x that we obtain 

from the Euler Square of index n, k. Without loss of 

generality, assume that 2^ < pI^ < pf^ < ... < p\\ 

let ki — max{2’',p{\p2^,... ,pp} = Clearly, 

ni = ^ = 2^p^f p'f ... p^'zf. Since the Euler Square 
of index ni,k exists and using it we construct nik x n1 
matrix, say with coherence at most -i. Now using these 
two matrices and a new matrix is obtained as 
detailed below: 

The matrix may be viewed as a block matrix possessing 
k blocks with each block being of size n x n^. Since ki > k 
and nifc = we have nik < n. The suitably zero-padded 
may be added as an additional column block to to 
generate a matrix of size nk x (n? + n{). One may obtain 
k number of zero-padded column blocks by placing in 
the corresponding row-block-locations of Using all these 
column blocks, we obtain a new binary matrix of size nk x 
{kn\), say Erom Lemma |2] the inner-product between 

any two columns in and is at most 1. Eurther the 
inner-product between any two columns one from and 
other from is at most 1, this is because the k nonzero 
entries of a column from fall only in one of the blocks 
of which has only one nonzero entry. 

The elements of the matrix may be generated as 

follows: for f = 1, 2,..., fc; (t — l)nf + 1 < j < 1 < 


i < (t — l)n, = 0, also for {t — l)nf + 1 < j < 

tn\, {t — l)n -f 1 < i < {t — l)n + nik, 


1 if (aj) 


= 

r 11 


-(l + (t-l)n.) 


- 1-1 


= i — {1 + {t — l)n) 
(mod m) 


0 otherwise. 


(5) 


and for {t — l)nf + 1 < j < tnf, (t — \)n + nik + 1 < i < 
nk, = 0, where (aj) is the fc—ad in the Euler Square 
of index ni, k. 


Similarly, using = max{2'^,pf,p^^,... 

712 = ^7 a CS matrix of size n2k x n\ (say <I>(^)) is 
obtained. The matrix may be viewed as a block matrix 
possessing k blocks with each block being of size ni x n\. 
As shown in Eig. 1, the matrix ikC) contains k such $C) 
matrices, therefore there are k^ blocks each of size ni x nf 
in ikC). Since ki > k2 > k and n2k = ^k, we have 
n 2 k < ni. The suitably zero-padded may be added 
in these blocks as additional column blocks to as 
shown in Eig. 1 and generate a new matrix, say of 

size nk x k^n\. Continuing this way till Euler Square of 
index ni,k, a matrix, say of size nk x k^nf is obtained 
at the stage. Einally, concatenating ..., 

results in a matrix with coherence at most t and column size 

(n^ -I- nfk + n^k'^ + n^k^ -f ... -f ufk!") > ^ . 

This methodology is explained diagrammatically in Eigure U 
for I = 2, that is, n = 2^pf^p^f. 



([)io] ^|(l) q;i2) 


Eig. 1: Block diagram depicting the extension of the column 
size of CS matrix constructed from the Euler Square. The way 
the blocks, obtained from several Euler Squares, are arranged 
ensures that the column size gets enlarged without affecting 
the coherence of initial matrix. 


The column extended matrix has RIP compliance as sum¬ 
marized by the following theorem. 

Theorem 6. The matrix rk*^^) ... satisfies the RIP 
with 5k' = for any k' < k 1. 

Remark 4: Combining the results in Sections IIII-BI and 
IIITCl it can be concluded that the constructed matrix attains 
column size m x c{mp)^ for some c G [1,2). In particular, 
when m is power of a single prime (that is, m = p\ I > 2), 
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c = 1, and in other cases c can be made to lie in (1, 2), as 
justified in Section UlI-CI 

IV. Comparison with existing constructions 

As discussed already in the Section III-Bl existing methods 
provide binary CS matrices for specific set of integers. The 
present work, however, provides a procedure for constructing 
binary matrices for a large class of row sizes. From Remark 
1, the maximum possible column size is between and 
2v?. Since A: > 2, it follows that (nk)^ > 4n^ and hence 
M < rn?. In the constructions provided by lIU, Q, llT4l . llT5l . 
M < m? holds, albeit for specific set of values of m. While 
in the construction of lfT3]l . M > m? holds. It is known ll25l 
that the Welch bound is not sharp in this case. Consequently 
the gain obtained in terms of increased column size does 
affect the coherence, which inturn restricts the sparsity of the 
solution to be recovered. In addition to providing general row 
size, the present construction is simple in the sense that it 
does not involve function evaluations like in and gives 
matrices with small density, which support algorithms with 
low computational complexity. For example, to generate an 
Euler square matrix of size p,p—l, it is only required to store 
two cyclic permutations of length p and p — 1 respectively 
and for index it is sufficient to store at most ^ 

permutations ll24l . 

To the best of our knowledge, the constructions possessing 
sparsity k' = y/m (that is, coherence p = -^) exist for non¬ 
binary matrices with row size m being prime or prime power 
El, a, 10, 0. While in our construction methodology, if 
we use the Euler Squares of index p,p — j for j = 1, 2, then 
the size of the matrix $ is p{p — j) x p^ and coherence p 
is Since the row size m of $ is p{p — j), we have 

// = ■ Hence, this matrix provides guarantees for signals 

of sparsity up to k' = [y/ml. This is true for any Euler 
Square of index — j, where p is a prime number, i > 1, 
and j = 1,2. Hence we construct the binary matrices that 
provide guarantees for signals of sparsity up to fc' = \j\/rn\ 
for different class of row sizes such as p*(p* — j), where p is 
prime i > 1, and j = 1,2. 

Eor an arbitrary binary matrix, if the inner-product between 
any two columns is at most 1, every column contains fixed 
i\/m) number of ones (that is coherence is at most -^) and 
row size is m, then the maximum possible column size M = 
0{m) as mentioned earlier. 

Using Euler squares, we have constructed matrices that 
provide guarantees for signals of sparsity up to fc' = \j\/rn\ 
with M = 0{m). Through an extension of our construction 
methodology, it is possible to generate ternary matrices, for 
which k' = ^ /m holds (that is, coherence is in the order 
of ^^ column size M = 0{m^), which is 
explained in the following: 

Suppose $ is a matrix constructed from the Euler square of 
index p®, p®— j for i > 1, j = 1 or 2 and H a Hadamard matrix 
of size (p® — j) or (p® — j) -f 1. To get a ternary CS matrix 
(i.e. matrix containing 0,±1 as elements) 'F, for each column 
of $, we replace each of its 1-valued entries with a distinct 


row of H and 0-valued entries with a zero row of same size. 
The size of the matrix T* becomes p®(p® — j) x. p^®(p® — j). In 
this construction, it is very easy to check that the c oherence 
of the matrix T* is which is in the order of y/ 

(Welch bound) with column size M = 0{mi). 

With a view towards comparing the numerical perfor¬ 
mance against the standard Gaussian (with entries drawn from 
^(0, ^)) and Bernoulli (with entries 
with probability i) random matrices, binary matrices of size 
55 X 121 and 230 x 529 are generated using the Euler Squares 
of indices (11,5) and (23,10). The OMP algorithm is used 
to solve the Iq minimization problem for a signal x. Let x 
denote the recovered solution. Erom the reconstruction, the 
Signal-to-noise ratio (SNR) of x is computed using 

SNR{x) = 10. logio { 

Vik- a;|| 2 / 

Eor recovery at each sparsity level, 1000 input A:—sparse sig¬ 
nals X (the nonzero indices are chosen uniformly randomly and 
the nonzero entries are drawn with ^ V(0,1)) are considered. 
The recovery is considered to be good if SNR{x) > IQQdB. 
Simulation results on these matrices indicate (shown in Eigure 
lU that the Euler Square matrix gives better performance than 
the Gaussian and Bernoulli random matrices for certain higher 
sparsity levels. While for other sparsity levels, these matrices 
give the same recovery performance. It is observed that the 
CS matrices constructed from the Euler Squares of indices p, s 
for a prime p and s = [|J top—1, give relatively superior 
performance than the other choices of s against their Gaussian 
and Bernoulli counterparts. 

A. Phase transition 

The phase transition diagrams depict the largest k (with 
fixed m and M) for faithfully recovering A:—sparse vectors 
via (i-norm minimization. In Eigure [3 the region above the 
curve is the one in which successful reconstruction is not 
possible, while the region below consists of points at which 
successful recovery is possible. Again, recovery is considered 
to be successful if SNR is greater than IQQdB. Given a set of 
points 5 = ^, we have generated phase transition by finding 
the largest sparsity k such that the successful recovery of 90 
percent is achieved by considering the average recovery over 
1000 iterations at each point by using the matrix which gives 
that particular point. Eigure |3 provides phase transition for 
Euler, Gaussian and Bernoulli random matrices for different 
values of 5 = ^ which were given by the matrices of size 
m X M, where m = 22,33,44,55,66,77,88,99,110 and 
M=121. It may be inferred from Eigure[3that the Euler Square 
based recovery is very competitive. 

B. Reconstruction of images 

The efficacy of Euler Square based matrix is demonstrated 
using image reconstruction from lower dimensional patches, 
where the patches are generated via the sensing matrices (as 
explained in more detail in the next section and Eigure |3l. 
Since the CS theory relies on the concept of sparsity, the key 
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Fig. 2: Comparison of the reconstruction performances of Eu¬ 
ler Square based, Bernoulli random and the Gaussian random 
matrices when the matrices are of size (a) 55 x 121 (top plot) 
and (b) 230 x 529 (bottom plot). These plots indicate that 
the Euler Square based matrix shows superior performance 
for some sparsity levels, while for other levels all matrices 
result in the same performance. The x and y axes in both 
plots refer respectively to the sparsity level and the success 
rate (in % terms). Eor matrices in (a) and (b), the coherences 
0.5226, 0.34 (for Gaussian matrices) 0.52, 0.28 (for Bernoulli 
matrices) and 0.2, 0.1 (for Euler based matrices). The values 
of coherence of Euler matrices are equal to their theoretical 
bound LL = 

“ m 



Eig. 4: Orginal image 



(d) (e) 


Eig. 5: Eor the original image in EigurelH the images in (b) and 
(d) are those reconstructed via the Euler Square based matrices 
with down-sampling factors 2.6 and 1.6 respectively. The 
images in (c) and (e) are those obtained via the corresponding 
Gaussian matrices. This figure states that Euler Square based 
CS matrices provide competitive reconstruction performance. 
The corresponding reconstruction errors are reported in Table 

m 



Eig. 3; Comparison of the reconstruction performances of 
Euler Square based, the Gaussian random and Bernoulli ran¬ 
dom matrices through phase transition. This plot indicates 
that the Euler Square based matrices provide wider recovery 
region than their Gaussian and Bernoulli counterparts. The 
X and y axes in the plot represent (m/M) and (k/M) 
respectively. This plot is generated for M = 121, m = 
22,33,44,55,66, 77,88,99 and 110. 


property needed of the image is enough sparsity in its original 
or some transform domain. 

As compressed sensing allows for the reconstruction of 
sparse vector x from its linear measurements ^x, we demon¬ 
strate reconstruction performance via a medical image. The 
reconstructions shown in Eigure |5] correspond to different 
down-sampling factors, viz 2.6 and 1.6. The associated re¬ 
construction errors in term of SNR are shown in Table U Here 
by down-sampling factor, we mean the ratio of original patch 
size to reduced patch size which is same as (where, mx M 
is the size of the matrix used for projecting data to the lower 
dimensional space). Erom Eigure |5] and Table H] it may be 
concluded that the Euler Square based matrices provide better 
reconstruction performance. 

V. Content based image retrieval (CBIR) using 
Euler Square matrices 

In this section, the usefulness of binary matrices so con¬ 
structed for the CBIR of medical databases is demonstrated. 
As more and more hospitals use picture archiving and com- 
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Down-sampling 
factor (N) 

Euler recovery 
SNR error 

Gaussian recovery 

SNR error (Average error) 

4 

13.36 

13.44 

2.6 

16.44 

15.03 

2 

19.63 

18.14 

1.6 

20.61 

19.74 


TABLE I; A comparative error analysis of reconstruction 
by Euler based and Gaussian matrices for different down- 
sampling factors (^) 4, 2.6, 2, 1.6. The average error over 
1000 iterations is reported for Gaussian matrices. 


munication systems (PACS), the medical imagery world wide 
is increasingly acquired, transferred and stored digitally ll20l . 
The increasing dependence on modern medical diagnostic 
techniques like radiology, histopathology and computerized 
tomography has led to an explosion in the number of medical 
images stored in hospitals. Digital image retrieval technique 
is crucial in the emerging field of medical image databases 
for clinical decision making process. It can retrieve images of 
similar nature (like same modality and disease) and charac¬ 
teristics. A typical CBIR system involves 2 steps, namely: 1) 
feature extraction or dimensionality reduction and 2) retrieval 
of relevant images through similarity metric. 

In fields like tele-medicines E], im the dimensionality 
reduction (DR) is needed for the analysis and classification 
of medical images, which is followed by the reconstruction 
at diagnostic center. Though conventional DR approaches 
ll^ project data to very low-dimensional spaces than the 
CS based method, in general, they fail to provide faithful 
reconstruction from reduced dimension. The DR based on 
Principal Component Analysis (PCA) has the potential to 
provide reconstruction as well. Nevertheless, it may be noted 
that PCA and Compressed Sensing become effective under 
different sets of conditions ED. Most importantly, the PCA 
based CBIR is very time-consuming as it is a data-driven 
approach. In ina, a dictionary learning (DL) based CBIR 
approach is presented, which learns dictionaries in Radon 
transform domain. The sparsity seeking DL approaches typi¬ 
cally exploit the framework of under-determined setting and 
hence work on some implicit assumptions on database. In 
CBIR applications, when data bases are not big enough, the 
sparsity seeking under-determined frame work may not be 
deployed efficiently. In addition, when labeled data are not 
used (as is the case with present work), one may not have 
enough members in a cluster, which prevents the applicability 
of Dictionary Learning ffSll . The present CBIR application 
involving the dimensionality reduction of database members 
is not prone to these problems. 

In general the CBIR approaches exploit the properties 
of database members for feature extraction lEl, ESI- For 
example, the CBIR method in uses the distribution of 
edges, and hence appears to work well on databases where the 
members have pronounced edges. The sparsity based methods, 
on the other-hand, exploit the presence of inherent sparsity 
present in most of the databases. In what follows, relevant 
features are extracted by projecting image content into lower 
dimensional space. 

The database members are divided into smaller 


patches {Ii,p\ I = l,2,...,iV and p = 1,2,...,M'} of equal 
size. Here, M' stands for the number of patches being carved 
out of each of database members. The vectorized versions of 
are then decomposed into the wavelet domain to generate 
sparse vectors, say ^ for each patch. It is to be emphasized 
here that one may use any transformation that sparsely repre¬ 
sents each patch. A down-sampled copy of 7;^ is generated 
via the binary sensing matrix T as 7"^ = TI'i^. If 7/^ is 
sparse enough, 7/^ (and consequently 7; p) can be recovered 
from the reduced vector I['^ using li norm minimization 
technique. Einally, the feature vector of database member 
is obtained as {7"p| p = 1,2,3,..., M'}, which is very small 
in size compared to {7/ p|p = 1,2, Given the query 

image Q, its feature vector {Qp\p = 1,2,..., M'} is generated 
similarly. Eor retrieving the images of database that are similar 
in content with the query image, cross-correlation is employed 
as a similarity metric (other criteria are also possible). The 
proposed CBIR method is shown in Eigure |6] 


Database Images 


Query Image 






Eig. 6: Block diagram of the proposed CBIR method. 

The performance of the image retrieval task is measured in 
terms of recall R = NJNm and precision P = / {N^ + Nf) 

where Nm is the total number of actual (or similar) images, 
Nc is the number of images detected correctly, and Nf is the 
number of false alarms. A good performance requires both 
recall and precision to be high, that is, close to unity. Recall is 
the portion of total relevant images retrieved whereas precision 
indicates the capability to retrieve relevant images. 

A. Experimental results 

Eor evaluating proposed CBIR, a database of 821 images 
present in the form of 12 classes (Eigure |7]i has been chosen. 
These classes of images containing skull, breast, chest and 
hand etc are taken from the popular IRMA databas^E Another 

2 WWW.irma-project.org 
















240 are considered as query images. 

From the database, = 1,2,3,821; p = 

1, 2, 3,64} and {//" |Z = i, 2, 3,821;p = 1, 2, 3,64} 
have been generated. are obtained using binary and 

Gaussian matrices respectively for the purpose of comparison 
of their performances in CBIR. 

The query image is compared with all database members 
in compressed domain using cross-correlation as similarity 
measure for retrieving top 10 similar images. The precision 
and recall of binary and Gaussian based feature vector with 10 
query images per class are shown in Table HIl In the simulation 
work, Haar wavelets have been used for sparsely representing 
each patch. 



nrapirara 




Classes 

10 query images per class 

Euler matrix 

Gaussian matrix 

Precision 

Recall 

Precision 

Recall 

Cl 

81 

47.8 

90 

58.1 

C2 

too 

97 

100 

85.5 

C3 

98 

89 

97 

81 

C4 

89 

58.7 

92 

75.2 

CS 

76 

48.9 

74 

45 

C6 

95 

69 

98 

71 

C7 

77 

39 

66 

36 

C8 

91 

61 

92 

45 

C9 

98 

82 

95 

79 

CIO 

100 

69 

100 

61 

Cll 

100 

85 

100 

82 

C12 

100 

46 

94 

37 

Average 

92 

67.4 

91.5 

63.4 


TABLE II: Performances (in % terms) of the proposed CBIR 
method based on the Gaussian and binary matrices. 


of the confusion matrix indicate if each class is classified 
correctly. The corresponding nonzero column entries indicate 
the presence of mis-classification by the CBIR scheme. From 
Table m it may be concluded that the binary based retrieval 
gives the performance comparable to that of its Gaussian 
counter part. The retrieval performances by both matrices may 
be improved by fine-tuning the associated parameters (like 
patch size) and by taking different sparsifying basis. Since the 
main focus of the present work is to construct CS matrices 
of general size, we do not go into the details of improvement 
of performance of CBIR and the reconstruction of database 
member from the reduced dimensional vectors any further. 







Fig. 8: Confusion matrix of Euler based binary matrix based 
CBMIR. 


VI. Concluding Remarks 


Fig. 7: Some of the images from each of 12 classes of database. 

The classification performances of different classes by bi¬ 
nary and Gaussian matrices are shown in Figures [8] and |9] in 
terms of a confusion matrix. The confusion matrix gives the 
accuracy of the classification results. The diagonal elements 


In this paper, Euler Squares of index (n, k) are used to 
construct binary sensing matrices of size nk x n^. Using these 
and with the help of prime factorization, matrices of general 
row size with asymptotically optimal column size and good 
coherence have been constructed. Further, using block wise 
extension, the column size of the constructed matrices has 
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Fig. 9: Confusion matrix of Gaussian matrix based CBMIR. 


been enhanced. Simulation results on two test cases show that 
the generated matrices capture the support of the unknown 
vector and give better performance than Gaussian matrices. It 
has also been demonstrated that the proposed binary sensing 
matrices exploit the pattern of inherent sparsity present in 
natural images and project data into lower dimensional spaces 
in such a way that the reduced vectors are useful for the 
purpose of CBIR. 
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